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Abstract 



^^ , The conformal equivalence of fourth-order gravity following from a 



non-linear Lagrangian L{R) to theories of other types is widely known, 



t^^ , here we report on a new conformal equivalence of these theories to 

■*««« . theories of the same type but with different Lagrangian. 

Q^i For a quantization of fourth-order theories one needs a Hamil- 

^ . tonian formulation of them. One of the possibilities to do so goes 

back to Ostrogradski in 1850. Here we present another possibility: 



A Hamiltonian H different from Ostrogradski's one is discussed for 
the Lagrangian L = L{q,q,q), where d'^L/d{q)'^ ^ 0. We add a suit- 
able divergence to L and insert a = q and b = q. Contrary to other 
approaches no constraint is needed because a = 6 is one of the canon- 
ical equations. Another canonical equation becomes equivalent to the 
fourth-order Euler-Lagrange equation of L. 

Finally, we discuss the stability properties of cosmological models 
within fourth-order gravity. 

PACS numbers: 04.50 Other theories of gravitation, 98.80 Cosmology, 
03.20 Classical mechanics of discrete systems 



*Extended version of a lecture read at the International School-Seminar "Problems of 
Theoretical Cosmology" , 1-7 September 1997, Ulyanovsk, Russia 



1 Introduction 

The purpose of this review is to present some of the points which have been 
discussed in recent years in connection with cosmological models following 
from fourth-order theories of gravity. 

Usually, one applies some version of a conformal equivalence theorem 
between these theories and Einstein's theory with additional scalar fields. 
But there exists also another type of conformal equivalence: In [1], we have 
shown that for non-linear L = L{K), G = dL/dR ^ the Lagrangians L 
and L{R) with L = 2R/G^ - 3L/G'^ gij = G^ gij and R = 3R/G^ - AL/G^ 
give conformally equivalent fourth-order field equations being dual to each 
other. The proof represents a new application of the fact that the operator 
□ — ^ is conformally invariant. 

The Ostrogradski approach [2] to find a Hamiltonian formulation for a 
higher-order theory is the most famous (see [2-9]) but not the only method. 
In set. 2 we present an alternative Hamiltonian formalism for fourth-order 
theories. It systematizes what has been sporadically done in the literature 
for special examples. 

Set. 3 deals with fourth-order gravity following from a non-linear La- 
grangian L{R). The instability of these theories from the point of view of 
the Cauchy problem is discussed. 

Set. 2 applies to general higher-order theories, set. 3 to gravity, and both 
are combined to fourth-order cosmology in set. 4. We discuss known exact 
solutions under the stability criteria mentioned before. 

In the final set. 5 we list some open problems for future research. 
The rest of this introduction shortly reviews papers on higher-order theo- 
ries. Eliezer and Woodard [2] and Jaen, Llosa, Molina and Vives [3] represent 
standard papers for the generalization of the Ostrogradski approach to non- 
local systems and to systems with constraints (see also [4-6]) applying Dirac's 
approach. 

Let the Lagrangian L be a function of the vector q^ and its first n temporal 
derivatives g^, Qa, ■ ■ ■ , q^^- The Hessian is 

^''^ = o (n)o (n) (1-1) 

and the non-vanishing of its determinant defines the regularity of L. In the 
following we do not write the subscript a; one can think of q as being a 
point particle in a (one- or higher-dimensional) space. In the Ostrogradski 
approach, Q = q is taken as additional position variable. This leads to an 
ambivalence of the procedure, because it is not trivial to see at which places 



q has to be replaced with Q, cf. [7]. We prevent this ambivalence in our 
alternative Hamiltonian, cf. set. 2, by putting Q = q. 

Ref. [8] discusses higher-order field theories. The problem is the lack of an 
energy bound, typically two kinds of oscillators with different signs of energy 
exist. Usually, one restricts the space of initial conditions to prevent negative 
energy solutions. The authors of ref. [8] redefine the energy analogous to 
the Timoshenko model, so one gets a positive mechanical energy inspite of 
an indefinite Ostrogradski Hamiltonian, they write: "An appealing aspect 
of this approach is the absence of any constraint." So it has this property 
in common with our approach set. 2, but it is otherwise a different one. 
Another standard procedure [3,8,9] for dealing with higher-order Lagran- 
gians is to consider them as a sequence in a parameter e, so one can break 
the Euler-Lagrange-equation into a sequence of second order ones. In [9] this 
is called "reduction of higher-order Lagrangians by a formal power series in 
an ordering parameter." [9] deals also with the Lie-Konigs theorem: a local 
Hamiltonian is always possible, and they consider some global questions. 

Let us show the famous counter-example [10]: it is an example of a 
second order system not following from a Lagrangian. 

X + y = y + y + ex = 

It follows from the Lagrangian 

L = -[f -y^ + e{xy - xy - i;^)] 

for e 7^ and has no Lagrangian otherwise. We mention this example to 
show that the following recipe need not to work always. Recipe for higher- 
order theories: "Write down the Euler-Lagrange equations, break them into 
a sequence of second order ones by introducing further coordinates. Find 
Lagrangians for these second order equations." 

A powerful method for dealing with a classical Lagrangian 

L = lg^Jq'q'-Viq) (1.2) 

is given in [11]. The Euler-Lagrange equation to Lagrangian (1.2) reads 

q' + r],q^' = -g^'Vk 
and is fulfilled for geodesies in the Jacobi-metric 

9ij = {E-V)gij 



Remark: For constant potentials V this is trivial, for non-constant potentials 
the constant E must be correctly chosen to get the result, for E = V it breaks 
down, of course. 

Stelle [12] cites Ostrogradski [2] but uses other methods to extract 
different spin modes for fourth-order gravity. In [13], a regular reduction 
of fourth-order gravity similar to the method with an ordering parameter 
mentioned above has been proposed as follows: In the Newtonian limit one 
has 

A$ + /3AA$ = AnGp, 

then one restricts to solutions which can be expanded into powers of the 
coupling parameter (3. Argument: If /? is a parameter, this is well justified, 
if it is a universal constant, then this restriction is less satisfying. Comment: 
This restriction excludes the usual Yukawa-like potential ^exp(— r/y^), so 
one may doubt whether this method gives the right solutions. Let us further 
mention ref. [14] for non-local gravitational Lagrangians like L = RO~^R-\-A 
in two dimensions and refs. [15-18] for the linearized i?^-theory. 

To facilitate the reading of set. 2, we pick up the example eq. (5) of [5]: 

L = [q^ + 4qq^ + 4q^]e^i (1.3) 

The equation of motion is [5, eq. (6)] 

2g(^) + 12gg(3) + 9{qf + ISq^q = (1.4) 

A good check of the validity of the formalism is the following: For a constant 
c > and g > 0, each solution of 

q = — 2g^ + c\Jq 

is also a solution of eq. (1.4). 

By adding a divergence to eq. (1.3) one gets L = (g)^e^'^. The alternative 
formalism requires to use q^ = q and q^ = q as new coordinates. So we get 

L = {q^f exp(3g^) (1.5) 

Eq. (1.5) represents the ultralocal Lagrangian mentioned in [5]. It is cor- 
rectly stated in [5], that the alternative formalism does not work for this 
version eq. (1.5) of the system. This clarifies that the addition of a diver- 
gence to a higher-order Lagrangian sometimes influences the applicability of 
the alternative Hamiltonian formalism. So one should add a "suitable" total 
derivative to the Lagangian. "Suitable" means, that the space of solutions 



is the same at both sides, and that the relation between the various coor- 
dinates is ensured without imposing any constraints. It turns out, that the 
Lagrangian L differing from L, eq. (1.3) by a divergence only 

L = -[f + 6qq^ + 2qq^%^'' (1.6) 

does the job. Of course, the variations of L, L, and L with respect to q 
all give the same equation of motion (1.4). But only in version (1.6) the 
alternative formalism (insertion of the equation Q = q and then apply the 
usual formulas of classical mechanics - to avoid ambiguities with the square- 
sign we have replaced q^ by q and g^ by Q) leads correctly to the Hamiltonian 
[5, eq. (7)]: 

H = --{pP-3P^Q)e-^'' + Q^e^'i (1.7) 

It essentially differs from the Ostrogradski approach because terms only linear 
in the momenta do not appear. The integrability condition Q = q and the 
equation of motion (1.4) both follow from the canonical equations of eq. 
(1.7); no constraint is necessary to get this. 

Wagoner [19] showed conformal relations between different types of sec- 
ond order theories (Einstein, Brans-Dicke), whereas refs. [20-23] discussed 
such relations between fourth-order and second order theories. Refs. [19-23] 
also discuss which of the conformally equivalent metrics can be considered 
physical. 



2 The alternative Hamiltonian formalism 

Let us consider the Lagrangian 

L = L{q,q,q) (2.1) 

for a point particle q(t), a dot denoting ^ and 

(„) ^ ^ 

The corresponding Euler-Lagrange equation reads 

Q^dL_ldL^^dL 
dq dt dq dt"^ dq 



We suppose this Lagrangian to be non-degenerated, i.e., L is non-linear in q. 
The highest-order term of eq. (2.2) is 



dW 



therefore, non-degeneracy (= regularity, cf. eq. (1.1)) is equivalent to require 
that eq. (2.2) is of fourth order, i.e. 



a(q) 



(If g is a vector consisting of m real components then this condition is to be 
written as Hessian determinant.) 



If we add the divergence ^G(g, q) to L, we do not alter the Euler- 
Lagrange equation (2.2). Furthermore, the expression J^G is linear in g, 
and so its addition to L does not influence the condition of non-degeneracy. 
The addition of such a divergence can therefore simply absorbed by a suitable 
redefinition of L. 

In the next two subsections we add a special and a more general diver- 
gence to get a Hamiltonian formulation different from Ostrogradski's one. In 
Kasper [4] a similar consideration has been made at the Lagrangian's level. 
Subsection 2.1 represents only a special case of subsection 2.2, but we write it 
down, because it has the advantage that the formulas can be given explicitly, 
and so the formalism becomes more transparent. 

2.1 A special divergence 

The addition of the following divergence is no more done by a redefinition of 
L 

LMv = j^[f{q)qql f{q) y^O (2.3) 

and we consider L = L + Ldiv The Euler-Lagrange equation is again eq. 
(2.2). Using 

we get 

L = L + nq)eq + f{q)m' + qq^'^] (2.4) 

which contains third derivatives of q. 
We introduce new coordinates 

a = q , b = q (2.5) 



(In the Ostrogradski approach, the second coordinate is q, instead.) It is 
obvious that there is exactly this one compatibility condition: 

'a = b (2.6) 

Let us insert eq. (2.5) into eq. (2.4). This insertion becomes unique by 
the additional requirement that L does not depend on second and higher 
derivatives of a and b, i.e., 

L = L{a,a,b,b) 

giving 

L = L{a, a, b) + f\a)d% + /(a) [b^ + db] (2.7) 

(In the Ostrogradski approach, there remains an ambivalence which of the 
q in the original Lagrangian is to be interpreted as second coordinate and 
which as time derivative of the first one.) 

The momenta are defined as in classical mechanics by 

dL dL 

(In the Ostrogradski approach, an additional term is necessary.) Inserting 
eq. (2.7) into eqs. (2.8) we get 

Pa = ^ + V\a)db + f\a)b (2.9) 

and 

Pb = /(a)d (2.10) 

Because of /(a) ^ 0, cf. eq. (2.3), we can invert eq. (2.10) to 

Pb 



Inserting eq. (2.11) into eq. (2.9) and dividing by /(a) we get 



(2.11) 



^ .[p„-^-2/'(a)6^] (2.12) 



/(a) da f{a) 

It is instructive to make a more general consideration: The question, whether 
eqs. (2.9, 10) can be inverted to a, b, can be answered by calculating the 
Jacobian (by the way: Carl Gustav Jacob Jacobi was born in Potsdam 10th 
December 1804.) 

J _ d{Pa,Pb) _ dpa dpb dpa Opb .^ ^^. 

d{a,b) da db db da 



We insert eqs. (2.9, 10) into eq. (2.13) and get 

J = - [f{a)Y (2.14) 

Because of / 7^ one has also J 7^ and the inversion is possible. This more 
general consideration gave the additional information that the Jacobian is 
always negative; this is one of the several possibilities to give the phrase 
"fourth-order theories are unstable" a concrete meaning. 
We define the Hamiltonian H as usual by 

H = apa + bph - L 

i.e., with eq. (2.7) we get 

H = dpa + bpb-L- f\a)a% - f{a)[b^ + db] (2.15) 

Here we insert a according to eq. (2.11) and get the Hamiltonian H = 
H{a,pa,b,pb). The factor of 6 in iJ automatically vanishes, so we do not 
need eq. (2.12). The canonical equations read 



dH 

dpa 

further 

dH 

dpb 

and 

dH 

da 

and 



d (2.16) 

b (2.17) 

-Pa (2.18) 



f)H 

The whole procedure is intended to give the following results: The Hamil- 
tonian H shall be considered to be a usual Hamiltonian for two interact- 
ing point particles a{t) and b{t). One of the canonical equations shall be 
equivalent to the compatibility condition eq. (2.6) and another one shall be 
equivalent to the original Euler-Lagrange equation (2.2), whereas the two 
remaining canonical equations are used to eliminate the momenta pa and pb 
from the system. The next step is to find those Lagrangians L which make 
this procedure work. From eqs. (2.15) and (2.11) we get 

rr VaVb j, Pb ,. vlf\a)b . .2 /o on^ 



In this form, eq. (2.16) coincides with eq. (2.11) and (2.17) with (2.12). So 
we may use eqs. (2.9, 10) in the following, because they are equivalent to 
eqs. (2.11, 12). 

Now, we use eqs. (2.19), cancel pf, by use of eq. (2.10) and get 

= ^ + 26/(a)-o/(a) (2.21) 

In order that the compatibility relation eq. (2.6) follows automatically from 

eq. (2.21), one has to ensure that /(a) 7^ (which is already assumed) and 
that 

identically takes place. The condition of non-degeneracy, 

is then also automatically fulfilled. One has the following possible Lagrangian 

L = --f{a)b'^ + K{a,a) (2.22) 

where K is an arbitrary function, but, for simplicity, we put K = 0. 

The last of the four canonical equations to be used is eq. (2.18) reading 
now with eqs. (2.9, 10, 20) 

= fb + 2/' db + -f b^ + f" a% (2.23) 

If we insert here eq. (2.5) we get exactly the same as the Euler-Lagrange 
equation (2.2) following from the Lagrangian 

L = -\f{qm' (2.24) 

Result: For every Lagrangian of type (2.1) which can be brought into 
type (2.24) with / 7^ the addition of the divergence (2.3) makes it possible 
to apply the new coordinates (2.5). Then the system becomes equivalent to 
a classical Hamiltonian of two particles, and the relation (2.6) between them 
follows without imposing an additional constraint. 



2.2 A general divergence 

In this subsection we try to generalize the result of the previous subsection 
by avoiding to prescribe the special structure (2.3) of the divergence to be 
added. We substitute eq. (2.3) by 

Ldiv = -^Hq, ^, q) (2.25) 

Keeping eqs. (2.5) we get instead of eq. (2.7) now 

L = L{a,d,b) + hid + h2b + hsb (2.26) 

where /i„ denotes the partial derivative of h with respect to its nth argument. 
Using eqs. (2.8), (2.10) is now replaced with 

Pb = hs{a,a,b) (2.27) 

Eq. (2.13) is kept, and (2.14) is replaced with 

J = -(/^23)' (2.28) 

We have to require that /i23 7^ 0, and then the equation pi, = h^ is locally 
invertible as d = F{pb,a,b). From this definition one immediately gets the 
identity Fi /i23 = 1. Two further identities to be used later are not so trivial 
to guess. To derive them, let us for a moment fix pi, and then calculate the 
increase of h^ with increasing a and b resp. The assumed constancy of /13 
yields the equations 

hi3 + F2h23 = (2.29) 

and 

h33 + F3h23 = (2.30) 

resp. to be used for deducing the generalization of eq. (2.21). One gets the 
result: For /i23 7^ (which is already presumed), the compatibility relation 
(2.6) follows automatically from the canonical equation (2.19) if and only if 

= L3 + h2 (2.31) 

is identically fulfilled. One can see: The condition of non-degeneracy of the 
Lagrangian (2.1) namely 

L33 7^0 

is equivalent to the condition /i23 7^ 0. For any given non-degenerate La- 
grangian we can find the appropriate divergence by solving eq. (2.31) as 
follows 

^(g,g,g) = -/ L3{q,x,q)dx (2.32) 

Jo 

10 



All other things are fully analogous: 

H =[p^- hi{a, F, b)]F - h2{a, F, b)b - L{a, F, b) (2.33) 

where F = F{pi,,a,b). Eq. (2.19) with (2.30) gives the compatibility con- 
dition (2.6). Eq. (2.18) with (2.29) is equivalent to the Euler-Lagrange 
equation (2.2). 

Let us summarize this section: For the Lagrangian L = L{q, q, q) where 
d'^L/d{qY 7^ we define L = L + Ldiv where 



d f dL 
Ldiv = --t: J -^{q,x,q)dx (2.34) 



We insert a = q and b = q, define the momenta Pa = ^ and P6 = ff and get 
the Hamiltonian H = apa + bpb — L. One of its canonical equations is a = 6 
and another one is equivalent to the fourth-order Euler-Lagrange equation 
following from L. By these properties, Ldiv is uniquely determined up to the 
integration constant. Contrary to other approaches, no constraint is needed. 



3 Instability of i?^-theories 

This section deals with the classical instability of fourth-order theories fol- 
lowing from a non-linear Lagrangian L{R). 

Teyssandier and Tourrenc [24], cf. also [25], solved the Cauchy-problem 
for this theory, let us shortly repeat the main ingredients. 

The Cauchy problem is well-posed (a property which is usually required 
to take place for a physically sensible theory) in each interval of -R-values 
where both dL/dR and d'^L/dR^ are different from zero. The constraint 
equations are similar as in General Relativity: the four Oi-component equa- 
tions. What is different are the necessary initial data to make the dynamics 
unique. More exactly: Besides the data of General Relativity one has to 
prescribe the values of R and ^ at the initial hypersurface. This coin- 
cides whith the general experience: Initial data have to be prescribed till the 
highest-but-one temporal derivative appearing in the field equation (here: 
fourth-order field equation, ^ contains third-order temporal derivatives of 
the metric). Under this point of view, classical stability of the field equation 
means that a small change of the Cauchy data implies also a small change of 
the solution. 

Now we are prepared to classify the stability claims found in refs. [26- 
34]. To simplify we specialize to the Lagrangian L = R — eR^ with the 
non-tachyonic sign e > and restrict to the range ^ > 0, i.e. R < ■^. 

11 



On the one hand, refs. [26,27,28] find a classical instability of the Minkowski 
space-time for this case of fourth-order gravity. Mazzitelli and Rodrigues [29] 
cite ref. [30] with the sentence "The Minkowski solution in general relativity 
has been proven to be stable." which refers to the positive energy theorem 
of general relativity. 

On the other hand, refs. [31-34] find out that the Minkowski space-time 
is not more unstable in this type of fourth-order gravity than in General 
Relativity itself. What looks like a contradiction from the first glance is 
only a notational ambivalence as can be seen now: The main argument in 
refs. [26-29] is that an arbitrarily large value ^ is compatible with small 
values of H^ and R^ . In refs. [31-34] however, following the Cauchy-data 
argument [24,25], (^ being part of the Cauchy data which are presumed 
to be small) stability of the Minkowski space-time is obtained in the ver- 
sion: If the Cauchy data are small (meaning: close to the Cauchy data of 
the Minkowski space-time) then the fourth-order field equation bounds the 
solution to remain close to the Minkowski space-time. 

The argument of ref. [31] is a little bit different: There the conformal 
transformation to Einstein's theory with a scalar field $ [20] is applied; it 
is observed that in the F(i?)-theory there are never ghosts which implies 
stability. Now, $ and ^ belong to the Cauchy data which is equivalent to 
the data i?, ^ in the conformal picture thus supporting the Cauchy data 
argument given at the beginning of this section. 



4 Cosmology 

Several papers [35-43] apply the conformal transformation theorem [20] to 
cosmology; so for interpreting the cosmological singularity [35], for deal- 
ing with anisotropic models [36,37], with transformation to Brans-Dicke ex- 
tended infiation [38]. Other papers apply this theorem as a mathematical 
device to transform exact solutions of one of the theories to solutions of the 
other theory. 

Chimento [44] found an exact solution for fourth-order gravity in a spa- 
tially fiat Friedmann model. He also found out that in the tachyonic-free 
case the asymptotic matter-dominated Friedmann solution is stable, and 
no fine-tuning of initial conditions is necessary to get the final (oscillating) 
Friedmann stage; particle production of non-conformal fields may backreact 
to damp the oscillations. [45] generalizes [44]: here the Dirac equation is 
considered, the result is that there appear also spinor field oscillations. 

Let us present the exact solution of [44]. For the spatially fiat Friedmann 
model with Hubble parameter H = a/a he solves the fourth-order field 
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equation with vacuum polarization term. The zero-zero component equation 
reads 

2HH -H^ + 6H^H + -H^ + H^ = (4.1) 

4 

The if ^-term stems from the vacuum polarization and the H'^-term. from the 
Einstein tensor. The remaining ingredients of eq. (4.1) come from the term 
R^ in the Lagrangian. (Here we only present the tachyonic-free case with 
A = and | in front of H^.) The factor in front of H^ should not influence 
the weak-field behaviour because for H ^ this factor only changes the 
effective gravitational constant. 

From eq. (4.1) the discussion of section 3 becomes obvious: (4.1) repre- 
sents a third-order equation for the cosmic scale factor a; it is a constraint 
and not a dynamical equation. (It is only due to the high symmetry, that 
accidentally the validity of the constraint implies the validity of the dynami- 
cal equation.) Supposed, eq. (4.1) would be the true dynamical equation for 
a theory, then the instability argument of [26-28] could apply. 

The ansatz for solving eq. (4.1) 

3s 
leads to a non-linear third-order equation for s 

2ss(^) -s^ + s^ = (4.2) 

Derivative with respect to t yields the equation s*-^-* + s = being linear in s 
and having the solution 

s = Ci + C2t + C3 sin(t + C4) 

Inserting this solution into the original equation gives the restriction |c2| = 
|c3|. Let us discuss this solution: C2 = leads to the uninteresting flat space- 
time. So, now let C2 7^ 0. Adding tt to C4 can be absorbed by a change of 
the sign of C3. Therefore, C2 = C3 without loss of generality. Multiplication 
of s by a constant factor does not change the geometry, so let C2 = 1. A 
suitable time-translation leads to ci = 0. Finally, the cosmic scale factor is 
calculated as a = s^/^ leading to 

2 
a = [t + sin(t + C4)]^/^ ~ t^/^ [1 + — sin(t + C4)] (4.3) 

The r.h.s. of eq. (4.3) gives the late-time behaviour deduced in [46]. The 
factor 1/t in front of the "sin" -term shows that the oscillations due to the 
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higher-order terms are damped. The total energy "sitting" in these oscilla- 
tions, however, remains constant in time (because of the volume-expansion), 
cf. [27], and can be converted into classical matter by particle creation. 

Let us mention some further cosmological solutions with higher-order 
gravity: [47] discusses the L(i?)-stability with a conformally coupled scalar 
field. Ref. [48] (partial results of it can be found in [49]) deals with fourth- 
order cosmological models of Bianchi-type I and power-law metrics, i.e. 

i=l 

with real parameters pi. The suitable notation 

3 

1=1 

gives the following: ai = a2 = 1 is the usual Kasner solution for Einstein's 
theory. a\ + a2 = 2ai is the condition to be fulfilled for a solution in L = 
R^. Refs. [50, 51] also discuss i?^-models. Ref. [33] considers inflationary 
cosmology with a Lagrangian 

L = R+ XR^^R^"'/R. 

Ref. [52] deals with anisotropic Bianchi-type IX solutions for L = R^. They 
look for chaotic behaviour analogous to the mixmaster model in Einstein's 
theory. Ref. [53] gives exact solutions for L = R^ and a closed Friedmann 
model, ref. [54] discusses the bounce in closed Friedmann models for L = 
R — eR^. Supplementing the discussion of [54, eq.(l)] let us mention: In 
the non-tachyonic case, there exist periodically oscillating models with an 
always positive scale factor a. Ref. [55] looks for chaos in isotropic models, 
e.g. by conformally coupled massive scalar fields in the closed universe. The 
papers [56,57] consider the stability of power-law inflation for L = R™- within 
the set of spatially flat Friedmann models. Refs. [58] give overviews on 
higher-order cosmology, especially chaotic inflation as an attractor solution 
in initial-condition space. [59] deals with quantum gravitational effects in the 
de Sitter space-time, and [60] gives a classification of inflationary Einstein- 
scalar-field-models via catastrophe theory. Ref. [61] considers Chern-Simon 
terms in Bianchi cosmologies and the cosmic no-hair conjecture. 



5 Summary 

The scope of this paper was to review recent work connected with higher- 
order theories, especially fourth-order gravity theories and their application 
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to cosmology. We presented a step necessary to deduce the Wheeler-de Witt 
equation for a cosmological minisuperspace model in fourth-order gravity and 
discussed the several stability claims and conformal transformation theorems. 

The method (first used in [51] for L = R^ and a spatially flat Friedmann 
model) to handle with eqs. (1.3 - 1.6) was systematically generalized in 
set. 2 (cf. also [62]) to give a Hamiltonian formulation of a general fourth- 
order theory. The possibility of deducing this method makes it clear that 
the method of ref. [51] is not restricted to highly symmetric models. The 
alternative Hamiltonian formulation has some advantages in comparison with 
Ostrogradski's one: No constraint is needed, the Hamiltonian is typically a 
quadratic function in the momenta. (Ostrogradski's approach leads always 
to a Hamiltonian linear in the momenta which gives artificial factors i in the 
Schrodinger equation.) The calculation of the momenta from the Lagrangian 
follows the usual equations (2.8) whereas the Ostrogradski approach needs 
some additional terms. Our approach is less ambiguous, cf. eq. (2.7). 

The fact that the Jacobian eq. (2.28) is always negative excludes the 
possibility to get a positive definite Jacobi metric in eq. (1.2). This is 
one of the many possibilities to say what is meant by the phrase "fourth- 
order theories are always unstable". The Jacobi metric plays the role of 
the conformally transformed superspace-metric used in quantum cosmology. 
And here the circle can be closed: In Einstein's theory (both for Lorentzian 
and Euclidean signature of the underlying manifold) the superspace-metric 
has Lorentzian signature and cannot be positive definite. So we get once 
more the result: Fourth-order gravity contains instabilities, but for the non- 
tachyonic case oi L = R — eR^ only those which it has in common with 
General Relativity. 

Let us finish by mentioning some open questions worth being attacked in 
the future: 

For other types of fourth- and higher-order theories the singularity be- 
haviour is not yet understood. 

To decide the quantum instability of the Minkowski or de Sitter space- 
times in fourth-order gravity one must solve the corresponding Wheeler-de 
Witt equations. In [29] they are deduced for the spatially fiat Friedmann 
model and the Lagrangian L = R — eR^. 

In [50] it is mentioned that a classical theory with higher derivatives has 
instabilities: "At the quantum level, the difference is even more dramatic. 
Noncommuting variables in the lower-derivative theory, such as position and 
velocities, become commuting in the higher-derivative theory." Remark of 
U. Kasper to this sentence: "The uncertainty relation is primarily between 
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positions and momenta. If the momentum is independent of the velocity 
then commuting position and velocity need not bother." 

Refs. [62-64] discuss the conformal transformations between fourth-order 
gravity to Einstein's theory with a scalar field. Renter [62] proposed to use 
the notion "Bicknell-theorem" to state that conformal equivalence. Dick [64] 
represents a valuable update of the discussion which of the two conformally 
equivalent metrics shall be considered physical. 

Because of the frequent appearance of singular points of the corresponding 
differential equations it is especially useful to have a growing set of exact 
solutions (see also [65]) of the non-linear gravity models. 
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